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ABSTRACT 

We use data from the Herschel- ATLAS to investigate the evolution of the far-infrared- 
radio correlation over the redshift range < z < 0.5. Using the total far-infrared lu- 
minosity of all > 5<7 sources in the Herschel- ATLAS Science Demonstration Field and 
cross-matching these data with radio data from the Faint Images of the Radio Sky 
at Twenty-Centimetres (FIRST) survey and the NRAO VLA Northern Sky Survey 
(NVSS), we obtain 104 radio counterparts to the Herschel sources. With these data 
we find no evidence for evolution in the far-infrared-radio correlation over the redshift 
range < z < 0.5, where the median value for the ratio between far-infrared and radio 
luminosity, gm, over this range is qiR = 2.40 ± 0.12 (and a mean of <?ir = 2.52 ± 0.03 
accounting for the lower limits), consistent with both the local value determined from 
IRAS and values derived from surveys targeting the high-redshift Universe. By com- 
paring the radio fluxes of our sample measured from both FIRST and NVSS we show 
that previous results suggesting an increase in the value of <7ir from high to low redshift 
may be the result of resolving out extended emission of the low-redshift sources with 
relatively high-resolution intcrfcrometric data, although AGN contamination could 
still play a significant role. 

We also find tentative evidence that the longer wavelength, cooler dust is heated by 
an evolved stellar population which does not trace the star-formation rate as closely as 
the shorter wavelength <250 /im emission or the radio emission, supporting suggestions 
based on detailed models of individual galaxies. 

Key words: galaxies: evolution - infrared: galaxies - radio continuum: galaxies 



1 INTRODUCTION 

The relation between far-infrared emission and the ra- 
dio e mission from galaxies at all redshifts is surprisingly 
tight Jde Jong et al.|[l985l : iHelou et al.lll985l : I Condon! [1992 : 
lGarrettH2002T ) and leads to the conclusion that both trace re- 
cent star-formation, in the local and distant Universe. The 
far-infrared emission is believed to arise from the thermal 
emission of dusty clouds surrounding regions of star forma- 
tion, whereas the radio emission arises from cosmic-ray elec- 
trons accelerated in supernova remnants of the dying stars, 
which emit synchrotron radiation. 

It is however unclear why there should be such a cor- 
relation between the thermal far-infrared emission and the 
non-thermal radio emission over such a wide range of galaxy 
types and masses, from starburst systems to more normal 
galaxies. As a result of this, many models resort to a rel- 
atively significant amount of fine tuning, such as assuming 
a much stronger magnetic field t han what is estimated via 
minimum energy arguments (e.g. iThompson et al. 1 120061 ). A 
full discussion of such arguments is beyo nd the scope of this 
pap er; however, we refer to the reader to lLacki et al.l (|2009r ) 
and lLacki fc Thompson! (|2009r ) who provide a detailed dis- 
cussion of the various physical interpretations of the far- 
infrared-radio correlation (FIRC). 

Observationally, recent work has concentrated on ex- 
ploring the FIRC as a function of redshift, mainly be- 
cause of the preponderance of deep Spitzer and radio 
data over relat i vely s mall < 10 square degree areas (e.g . 
Appleton et al l 12004 iFraver et ail |2006|; llbar et all l20Qg|: 



Murphy et~all 12009: Michalo wski et al.l [20T0I : ISargent et all 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and 
with important participation from NASA 
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l2010l : iBourne et al.llioioh . This has led to several authors 
suggesting that the F IRC remains uncha nged out to high 
redshift (z>1.5) (e.g. ISargent et all l201(il ). whereas others 
suggest a shallow decrease in the rati o between far-infrare d 
luminosity and radio luminosity (e.g. ISevmour et ail feoog). 
Constraining the evolution of the FIRC is important as 
it may help our understanding of the physical mechanism 
which results in such a tight correlation between the ther- 
mal and non-thermal emission. For example, metallicity and 
temperature evolution would effect the far-infrared luminos- 
ity but not necessarily the radio synchrotron emission, where 
as an evolution in the magnetic field strength would alter 
the radio emission. Furthermore, evolution of the gas den- 
sity and the scale height of a galacti c disk may also influ- 
ence both emission mechanisms (e.g. IHelou fc Bicavlll993l : 
iNiklas fc Becklll997l ). 

The Balloo n-borne Large Aper ture Submillimetre Tele- 
scope (BLAST; iDevlin et alH iooa 1 ) experime nt followed by 
the la unch of the Herschel Space Observatory jPilbratt et al.l 
2010) open up the possibility of fully investigating the 
FIRC using spectral energy distributions which have been 
measured with far-infrared and sub-millimetre photome- 
try. The first of these concentrated on the relatively deep 
fields of the extended Ch andra Deep Field South with 
BLAS T Jlvison et al.ll2010al l and GOODS-North with Her- 
schel jlvison et al. 2010bl h Both of these studies find evi- 
dence for a modest evolution in the FIRC from low-redshift 
through to the high-redshift Universe, although if confined 
solely to the far-infrared selected galaxies in these samples 
the evidence for evolution is not as significant. 

However, we are still lacking a statistically significant 
sample of objects which can be used to determine the form 
of the FIRC in the low - redsh ift (0.05 < z < 0.2) Universe. 
lYun, Reddv fc Condonl (|200ll ) used data from the In frared 
Astronomical Satellite (IRAS; iNeugebauer et aT]|l984l ') com- 
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bined with the NR AO VLA Northern Sky Survey (NVSS; 
ICondon et al.| [l998) to investigate the FIRC in the very low- 
redshift Universe (i.e. z < 0.05) and found that in their 
sample of 1809 galaxies 98 per cent were consistent with 
a simple linear FIRC. The next step is obviously to push be- 
yond these very low redshifts and to fill the gap in the mea- 
surem ents of the FIRC between the lYun. Reddv fc Condor] 
l|200ll ) sample at z < 0.05 and those conducted on smaller 
fields at higher redshifts (z > 0.5). 

We are now in a position to do this with the data 
acquired as part of the Hersc hel Astrophysical Tera-Hertz 
Large Area Survey (H-ATLAS: lEales et all2010l ). H- ATLAS 
is the largest open time key project on the Herschel Space 
Observatory. The survey aims to map 550 deg 2 at five wave- 
lengths from 100 — 500 fim to 5<r point source sensitivities, 
including confusion noise, of 132, 126, 32, 36 and 45 mjy at 
100, 160, 250, 350 and 500 /im respectively. The main scien- 
tific goals are to gain a complete understa nding of galaxies 
in the relatively low -redshift Universe (e.g. iDve et afll2010l ; 
lAmblard et aljboioh . to study the evol ution of rare objects , 



such as active galactic nuclei (AGN; Serjeant et al.l |2010| ; 



Hardcastle et al.ll20ld ; lBonfield et alj|20ich and galaxy clus- 
ters, and also to act as a unique dat a set for finding distant 
lensed galaxies l|Negrello et alj|2010h . In this paper we mea- 
sure the FIRC using the Science Demonstration Phase of the 
H- ATLAS survey which covers ~ 14 degree 2 over an equato- 
rial field centred at 09 h 05 m 30 s , +00°30'00", corresponding 
to one of the fields within t he Galaxy and Mass Assembly 
(GAMA; iDriver et all 120091 ') Survey which has a wealth of 
other multi-wavelength data. 

The paper is arranged as follows: In Section [2] we dis- 
cuss how our sample was selected using far-infrared data 
from Herschel, along with optical, near-infrared and radio 
data. In Section we measure the FIRC and investigate 
how this may evolve with redshift and whether there is any 
dependence on either radio or far-infrared luminosity. In Sec- 
tion0]we present our conclusions. Throughout the paper we 
use a concordance cosmology with Ho = 70 km s _1 Mpc~ , 
Clyi = 0.3 and £2a = 0.7. 
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Figure 1. Integrated far-infrared luminosity versus redshift for 
the > 5ct sources from the H-ATLAS SDP field. The small dots 
denote all of the H-ATLAS sources which we use for the stacking 
of the radio data (see Section 12.411 , while the filled circles denote 
those sources with > 5<r radio detections for which we can calcu- 
late the FIRC directly. The open circles represent those sources 
with Li.4GHz > 10 23 W Hz -1 , i.e. in the AGN-luminosity regime, 
and the open stars are those objects with morphological signs of 
AGN activity in the radio maps (see Section [3). 



full catalogue - to determine the ratio of the probability that 
an optical source is related to the 250 /jm source to the prob- 
ability that the optical source is unrelated. One can then use 
these values to determine a reliability, R, which estimates 
the probability that a given source is the true counterpart 
to the Herschel object. The likelihood ratio accounts for the 
fact that not all Herschel sources will be detected in the 
SDSS r-band data, due to many galaxies residing at higher 
redshift than the SDSS is sensitive to, and also possibly due 
to dust obscuration. In what follows, we only consider those 
2334 sources with reliable counterparts (R > 0.8) to the 
//ersc/ieZ-selected catalogue. 



2 SAMPLE SELECTION 



For this study we use the H-ATLAS Science Demonstra- 
tion Phase (SDP) data. Full details o f the Spectral and 
Photometric Imaging Receiver (SPIRE; iGriffin et ai]|2010fi 



and Photodetector Array Camera (PACS; Poglitsch et al.l 



\201(j) data reduction for the H-ATLAS data can be found 



Pascale et all l|2010h and llbar et al.l (l2010li respect ively. 



The H-ATLAS 5a source catalogue (|Rigbv et alj|2010h was 
subsequently cross- matched with the Sloan Digital Sky Sur- 
vey Data Release 7 (SPSS; IXbazaiian et alj 120091 ) and the 
United Kingdom Inf rared Deep Sky Survey Large Area Sur- 
vey (UKIDSS-LAS; lLawrence et all 120071 ). F ull details of 
the cr oss-matching procedure can be found in ISmith et al.l 
(2010); however to summarize, we use the SDSS DR7 r-band 
catalogu e, limited to r ^ 22.4, and perf orm a likelihood ratio 
analysis (jSutherland fe Saundera ll992') on all possible coun- 
terparts within 10 arcseconds of the 250 /J,m positions (the 
beamsize at 250 /im is 18 arcseconds). The likelihood ratio 
method relies on using not only positional offsets, but also 
magnitude information - both of the counterparts and of the 



2.1 Spectroscopic and Photometric redshifts 

For our study on the evolution of the FIRC, redshift in- 
formation for all of the sources is crucial. For the brighter 
objects in our sample we are able to use the SDSS spectro- 
scopic redshifts, for the fainter sources we use the year 1 
and 2 d ata from the Gala xy and Mass Assembly (GAMA) 
survey |Priver et alj|20ioh where available. The target se- 
lection for the GAMA survey was r < 19.4 or z < 18.2 
or Kab < 17.6, full details of th e target select i on an d the 
target priorities can b e found in Bal drv et al.l (|201dh and 
iRobotham et all (|2010| ) respectively. The combination of the 
SDSS and GAMA redshift surveys results in 877 sources 
having spectroscopic redshifts. 

Using the combination of SDSS ugriz and UKIDSS- 
LAS YJHK data we determine photometric redshifts 
for the sources without spectro scopic redshifts. We used 
the publicly available ANNZ l|Collister fc Lahavl l2004f ) 
photometric redshift code, using spectr o scopic redshifts 
from the SPSS, GAMA (|Priver et all l2010l), D EEP2 
dDavis et all l2007h, zCOSM OS (|Lillv et al.l l2007h and 
2SLAQ ( Cannon et alj |2006). This allows us to construct 
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a spectroscopic training set with large numbers of objects 
(> 1000 per bin of unit magnitude or 0.1 in redshift) up to 
r-band magnitudes r < 23 and redshifts z < 1.0, i.e. to ap- 
proximately the photometric depth of SDSS and UKIDSS- 
LAS. Fu rther details of the photometric redshifts can be 
found in ISmith et all l|2010h . For our analysis we always 
use the spectroscopic redshift, where available, in prefer- 
ence to the photometric redshifts which have a typical error 
of Az/(1 + z) ~ 0.03. In all cases where we has a spec- 
troscopic redshift, this agreed with the photometric redshift 
within the uncertainties. 



2.2 The radio data 

Based on this Herschel catalogue, with optical identifica- 
tions, we measured radio fluxes directly from the imaging 
data of both the NVSS and the Faint Imag es of the Ra- 
dio S ky at Twenty-one centimetres (FIRST; iBecker et al.l 
Il995h at the position of the optical galaxy. The NVSS has 
a typical rms flux-density of 0.45 mjy and has the benefit 
that it was carried out in VLA D-Array and therefore has 
a resolution of around ~ 45 arcsec which is larger than the 
expected size of any individual galaxy within the H- ATLAS 
galaxy sample. As such we expect little flux to be missed or 
resolved out using the NVSS. On the other hand the FIRST 
survey with a 5-arcsec beam is deeper than NVSS (rms = 
0.15 mjy) but may miss extended flux, particularly at the 
low-redshift end of our source distribution. We therefore use 
the NVSS survey to measure the radio fluxes at the positions 
of the H- ATLAS sources with spectroscopic or photometric 
redshifts z < 0.2 and the FIRST survey to measure the ra- 
dio fluxes of sources with z > 0.2, where the FIRST beam 
equates to a galaxy size of >16 kpc. We include all sources 
with a flux-density 5i.4GHz > 5 times the background noise 
level determined using each individual 5x5 arcmin 2 cut-out 
around each source. We note that only two of the sources 
with direct radio detections at > 5a lie above a redshift of 
z — 0.2 Therefore this split makes little difference to the 
results presented in this paper, but may be important when 
comparing to results which investigate the FIRC at higher 
redshifts. To summarise, our sample is comprised of all ob- 
jects detected at ^ 5a in the H- ATLAS catalogue, with radio 
counterparts detected at > 5a above the local background 
rms from the NVSS (z < 0.2) and FIRST (z > 0.2) ra- 
dio surveys. It is possible that the lower resolution NVSS 
survey could suffer from some source confusion, where mul- 
tiple sources contribute to the measured flux of the object. 
We checked for this using the higher resolution FIRST data 
to identify whether there was more than one source within 
the NVSS beam, and where this was deemed to be an issue 
(< 5 per cent) the FIRST maps was used to measure the 
flux. We also checked for extended emission associated with 
AGN as detailed in Section 



2.3 Far-infrared luminosities 

Far-infrared luminosities were derived using two methods. 
The first uses a simple modified blackbody fit to the three 
SPIRE bands at 250, 350 and 500/im, and the PACS data at 
100 and 160/^m where available. We fixed the dust emissivity 
index to f3 = 1.5 and varied the temperature over the range 



10 < T < 50 K, following iDve et all (|2010l ). As in lDve et al.l 
(2010) we find a median temperature of T = 26 K. For each 
source we measured the integrated far-infrared luminosity 
(8 — 1000 (im) and the rest-frame 250 /im luminosity. We 
also fit the data with an emissivity index of /3 = 2, which 
gave a median tempe rature of T = 23 K, again consistent 
with IDve et al l l|2010l ). However, the integrated far-infrared 
luminosities derived for both fits were consistent within the 
uncertainties with a median offset of +0.07 dex in integrated 
far-infrared luminosity for /3 — 1.5 compared to those fit 
with /3 = 2. 

It is widely acknowledged that simple modified black- 
body fits may underestimate the total far-infrared emission, 
particularly at the short wavelength end of the spectrum. 
For this reason we also determine the far-infrared lumi- 
nosity of our so urces using the comprehensive set of star- 
burst SE Ds from Sicbcnmorg en fc Kriigell ([2007), in a simi- 
lar way to lSvmeonidis et al.l |2008l) . All star-burst and Ultra- 
luminous Infrared Galaxy (ULIRG) SEDs were considered 
to fit the datcQ- We use all three SPIRE data points and the 
PACS detections where available. In the absence of a detec- 
tion at the PACS wavebands we use the 2a upper limit to 
constrain the SED to lie below these values. We note that 
no shorter wavelength data is available for the H-ATLAS 
field at the present time. We find that the templates pro- 
vide far-infrared luminosities approximately 0.7 dex more 
luminous than the corresponding modified blackbody fits, 
although the monochromatic luminosities at the rest-frame 
wavelengths covered by the SPIRE data are consistent, as 
expected. For the purposes of this paper we use the far- 
infrared luminositie s determined from the ful l SED fits using 
the templates from ISiebenmorgen fc Kriigell (|2007l ) and for 
the monochromatic luminosities we use the modified black- 
body fits, although we emphasise that these are consistent 
between the two methods. 

Fig. [T] shows the distribution of our sources on the 
far-infrared luminosity versus redshift plane, which demon- 
strates that we are able to sample the z < 0.2 Universe suffi- 
ciently well for individual sources to investigate the form of 
the FIRC with this sample. We use stacking of the radio data 
(Sec. 12. 4p to determine whether our findings at z < 0.2 are 
consistent with the average sources properties up to z ~ 0.5 
and beyond, by comparing with previous results in t he lit- 
erature (e.g. Ilvison et al.ll2010aH bl; iBourne et al.ll201(il). 

W e also use the results described in lHardcastle et all 
(2010) to obtain the far-infrared luminosities for those 
sources which are strongly detected at radio wavelengths 
but which fall below the 5a detection limit of our Herschel 
catalogue, and thus do not make it into our sample. The 
combination of these allows us to probe to lower and higher 
radio luminosities at all redshifts, thus negating some of the 
biases inherent to selecting in a single band. We show these 
for completeness; however, we do not use them in calculating 
the FIRC. 

In Fig. [2] we show the correlation between the far- 
infrared luminosity and the rest-frame 1.4 GHz radio lu- 
minosity. We assume a radio spectral index of a = O.£0, as 



1 see http://www.eso.org/~rsiebenm/sb_models 

2 Wc use the definition of spectral index such that the radio flux 
density S v oc v~ a , where v is the frequency. 
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Figure 2. Integrated far-infrared luminosity versus rest-frame 
1.4 GHz radio luminosity of the high-reliability 250 (im detected 
sources in the H- ATLAS SD field with radio cross-identifications 
> 5ct above the local noise in the NVSS and FIRST images. NVSS 
fluxes were used for the sources at redshifts z < 0.2 and FIRST 
for those sources with z > 0.2, as explained in the text. The open 
stars represent those sources with morphological signs of AGN 
activity (see Section [3]l and the open circles denote those sources 
with Li.4GHz > 10 23 W Hz -1 , and thus have a high probability 
of AGN contamination in their radio emission. The arrows denote 
the 5<t upper limits in the radio luminosity from the NVSS and 
FIRST. 



found for submillimctre galaxies by llbar et alJ fcoid ). for 
the 104 objects (72 with spectroscopic redshifts and 32 with 
photometric redshifts) which lie above our flux limit in H- 
ATLAS and with > 5 a detections in the radio maps. One 
can easily see the strong correlation between the far-infrared 
and radio emission in our sample, with the best fit relation 
log 10 Lpir = (17.8 ± 0.7) + (0.42 ± 0.07) log 10 Li. 4G Hz for 
those sources with L\a ghz < 10 23 W Hz" 1 (see Section EJ). 



2.4 Stacking 

For those sources detected at Herschel wavelengths but 
which fall below our 5a threshold from the radio data, we 
also perform a stacking analysis on the NVSS and FIRST 
radio images to obtain sta tistical detec t ions. To do this we 
followed the technique of IWhite et al.l l|2007h . which clips 
the image at 5a around the median value and then com- 
bin es the individual im ages weighted by their variance (see 
also lFalder et al . 2010). The stacked images are generated in 
various bins of far-infrared luminosity and redshift to allow 
direct comparisons to those sources with direct detections. 
Again we do not use these to derive the median and mean 
values for the FIRC but just as a guide so the reader can 
judge for consistency. 



3 THE FAR-INFRARED RADIO 
CORRELATION 

The FIRC we use here is defined as the logarithmic ra- 
tio of the integrated far-infrared flux, Sir, deter mined from 
rest-frame wavelengths 8 — 1000 fim (following Belli 120031 ; 



llvison et all2010al ). and the rest-frame 1.4 GHz fc-corrected 
flux density, Si.4GHz such that gm, = log 10 [(5m/3.75 x 
z )] and is dimensionless. 5ir has units of 
W m -2 , 3.75 x 10 12 Hz is the normalising frequency (Helou 
et al. 1985) and 5i,4 GHz has units of W m~ 2 Hz -1 . Following 
llvison et all ri2010ah . we also calculate the monochromatic 
relationship between the fc-corrected emission at 250^im (us- 
ing our modified blackbody fit) and 1.4 GHz (using a spec- 
tral index of a = 0.8), hereafter g250- 

Contamination in the sample from radio-luminous AGN 
would preferentially decrease the value of giR. The ra- 
dio luminosity at which AGN begin to dominate the 
source population occurs around Li.4GHz~1 23 WHz -1 (e.g. 
iMauch fc Sadlerll2007l ; IWilman et all 120081 '). Inspecting the 
FIRST radio images we find that nine of the sources in our 
sample show clear signs of AGN activity, i.e. jet-like phenom- 
ena. All of these sources have radio luminosities Li.4ghz > 
10 23 W Hz" 1 and qm. values below the sample mean, as 
would be expected for an AGN dominated source. We omit 
these nine sources from all subsequent analyses, even though 
the effect they would have on our results is negligible. All 



objects with radio luminosities L1.4 ghz > 10 
represented with open symbols in all figures. 



WHz- 1 are 



3.1 Redshift dependence of the far-infrared— radio 
correlation 

In Fig. [3]we show the giR parameter as a function of redshift 
for all of the objects in our sample with > 5a detections at 
1.4 GHz. If we consider this entire sample then we obtain a 
mean value of giR = 2.15±0.09 and median giR. = 2.22±0.11. 
We also measure the mean and median values of giR after ex- 
cluding the radio sources with obvious signs of AGN activity 
and those above L1.4 ghz > 10 23 W Hz -1 . Unsurprisingly 
this leads to a larger value of the mean, giR = 2.33±0.09 and 
median giR = 2.40 ±0.12. To estimate the effect of the lower 
limits in our measureme nt of the value for giR w e use the 
Kaplan-Meier estimator (|Feigelson fc Nelsonl fl985). which is 
a maximum likelihood estimator to calcu late the mean in th e 
presence of lower or upp limits (see e.g. lSaiina et alj|2008l ). 
Using this method we find a mean of giR = 2.52±0.03. These 
are consistent with the value of giR measure d in the very 
local Universe (|Yun. Reddv fc Condor] 120011 ) using IRAS 
data, who found giR = 2.34 ± 0.01, although we note that 
lYun. Reddv fc Condon! (2001) calculated the value of giR 
based on the 40 < A < 120 /im wavelength range only. 
iBourne et ail (|2010h also discuss this point and suggest that 
the difference between the value of giR measured in this lim- 
ited wavelength ranged compared to the total infrared lumi- 
nosity is ~ 0.3 dex for an M51 template SED, but note that 
this is subject to large variations due to different relative 
contribution from cooler dust. We use our models to find 
what range this correction spans and we also find a median 
difference of 0.3 dex with a standard deviation of 0.4 dex 
in the fra-infrared luminosity; thus the consistency remains 
at the 2 a level with the results of lYun. Reddv fc Condon! 
(2001). Comparing with the results of lBelll (120031 ) . who used 
the bolometric far-infrared luminosity of local galaxies, and 
found giR = 2.64, we find that our results are slightly lower, 
but again within 2 a. Our value of giR is also in good agree- 
ment with th e values measured a t higher redshift by va rious 
authors (e.g. llvison et all feoiOa b; Sargent et al.ll2010l ) who 
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Figure 3. The FIRC versus redshift for our complete sample. 
Solid circles denote the sources with Li.4GHz ^ 10 23 W Hz -1 , 
whereas the open circles represent the sources with £<i.4GHz > 
10 23 W Hz -1 , i.e. the radio luminosity at which AGN begin to 
dominate the source population. The dashed line is the weighted 
mean value of qm for the Li.4GHz < 10 23 W Hz -1 sources and 
the grey shaded region shows the standard deviation on qm- The 
arrows denote the 5<r lower limits based on the upper limits from 
the radio data. 

find values in the range 2.4 < qm < 2.8. Some find evidence 
for an increase in qm towards low redshift fe.g. lBourne et all 
2010), but our data do not support this; we return to this 
point in Section [IOTT] 

Fig. [4] shows the value of 5250 against redshift for our 
sample. Only considering those objects with Li.4GHz < 
10 23 W Hz -1 , we determine a mean 5250 = 1-78 ± 0.04 
an d median 325 = 1.8 4 ± 0.05, similar to the values quoted 
bv llvison et al.l (|2010al ) with the fe-correction applied. Using 
the lower limits and the Kaplan-Meier estimator we find a 
mean 5250 = 2.01 ± 0.04. However, we find that the disper- 
sion on (7250 is a factor of two tighter than the corresponding 
dispersion for qm- This suggests that the flux measured at 
rest-frame 250 /im is possibly more closely related to the 
radio emission than the integrated far-infrared luminosity, 
although we note that the larger dispersion could be due to 
the large uncertainties of the complete SED due to our lack 
of direct detections in the PACS observations, particularly 
for the sources lying close to the flux-density limit at SPIRE 
wavelengths. 

However, there are reasons why we might expect to 
see a larger dispers ion in qm. than in g 2 so • For example, 
iDunne et alj ((2000b) showed that the relationship between 
the radio emission and the 850/im emission measured with 
SCUBA has significantly more scatter than in the similar re- 
lationship between radio emission and the far-infrared emis- 
sion determined from a fit to the SED which includes data 
from IRAS. Similarly, IVlahakis et"afl (|2007l l showed that 
the relation between 850/im emission and radio emission 
has significantly more scatter than the relationship between 
60/im and 1.4 GHz. Both suggest that this could be due 
to the presence of cold dust which dominates the emission 
at sub-millimetre wavelengths, but which is not as strongly 
associated with massive star formation. Given that our to- 
tal far-infrared luminosities are largely constrained by the 



Figure 4. The g250-radio correlation as a function of red- 
shift. The open circles denote those sources with L1.4 GHz > 
10 23 W Hz -1 . The dashed line is the weighted mean value of 9250 
for the L1.4GHZ < 10 23 W Hz -1 sources, and the grey shaded re- 
gion shows the standard deviation on 9250 ■ The arrows denote the 
5cr lower limits based on the upper limits from the radio data. 

longer wavelength SPIRE data points, then we would ex- 
pect the scatter in these value to be larger than the scatter 
in monochromatic luminosities derived at 250 fim which in 
this scenario would more closely trace the dust heated by 
young stars. 

This explanation can also be reconciled with detailed 
studies of individual g alaxie s. Using Herschel observations 
of M81, iBendo et alj (|2010h recently found evidence that 
the longer wavelength (160-500/im) emission traces dust 
with temperature 15-30 K which is heated by evolved stars, 
whereas the shorter wavelength (~ 70/im) emission is domi- 
nated by dust heated by young stars in star- forming regions. 
Further more, in a sam p le of nearby galaxies observed with 
Spitzer, ICalzetti et alj (|2010l ) also found evidence that as 
one moves to longer wavelengths, the conversion to a star- 
formation rate becomes more uncertain with the dispersion 
increasin g by a factor of ~ 2 m oving from 24pm to 70^im, 
although lAppleton et al.l |2004h find the opposite to be true, 
this could possibly due to the wide redshift range, < 2 < 1 
(and thus luminosity range) in the Appleton et al. study. 
Thus, it appears that we may be seeing this effect on the 
FIRC, where the radio emission remains a robust tracer of 
the star-formation rate but the longer wavelength Herschel 
data has an increasing contribution to the emission from 
dust heated by evolved stellar populations , compared to the 
heating by you ng massive stars (see also iRowan-Robinsonl 
l2003l ; lBe3l2003h . 

3.2 Luminosity dependence of the 
far-infrared— radio correlation 

With the large number of objects available to us due to the 
high-sensitivity of the Herschel Space Observatory, we are 
able to investigate the form of the far-infrared radio correla- 
tion over a range of redshifts, far-infrared luminosities and 
radio luminosities. However, both the far-infrared data and 
the radio data are derived from flux-density limited par- 
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ent images and any analysis is susceptible to Malmquist 
bias, where the higher redshift sources are biased towards 
the most luminous sources, and Eddington-type bias for 
signal-to-noise limite d cat alogues. The r efore, f ollowing both 
ISareent et ail (|2010h and llvison et all l|2010bh we split our 
sample into three radio luminosity bins (Fig. [5| and three 
far-infrared luminosity bins (Fig. [6} and determine the mean 
and median value for giR, in each of these bins. 

3.2.1 Binning in radio luminosity 

Fig. [5] shows that there is very little evidence for strong 
evolution in the FIRC over the redshift range < z < 0.5 in 
any of the individual radio luminosity bins, where a constant 
value of giR is consistent with the data. We also find similar 
results when using 5250- However, we find that the mean 
(and median) values of gm (5250) decrease with increasing 
radio luminosity from gm = 2.55 ± 0.11 (5250 = 2.23 ± 0.10) 
for L1.4GHZ < 10 22 W Hz" 1 to g m = 2.35 ± 0.11 (g 25 o = 
1.93 ±0.04) for 10 22 < (L 1AGHz / W H Z - 1 )< 10 23 , through 
to qin = 2.33 ± 0.24 (g 250 = 1.81 ± 0.06) for Li. 4G Hz > 
10 23 W Hz -1 . The most straightforward explanation of this 
is that there is an increasing fraction of AGN as we push to 
higher radio luminosities. 

Another possibility is that the increase in gm and 5250 
towards low redshift could be caused by the increasing prob- 
ability of resolving out extended structure with relatively 
high-resolution radio data. We do not believe that this is an 
issue with our sample as we have restricted ourselves to using 
the low-resolution NVSS at z < 0.2 and only use the FIRST 
survey at the higher redshifts 2 > 0.2. However, this could 
offer an explanation of the increase in girt towards low red- 
shift seen with surveys targeting the high-redshift Universe. 
At high redshift (z > 0.5) short baseline observations are 
not as important, due to the fact that distant star-for ming 
galaxies are relatively compact (e.g. lMuxlow et al J2005h and 
are fully resolved with longer baselines, whereas short base- 
lines are needed for studying lower redshift galaxies which 
have larger angular extent on the sky. This extended emis- 
sion would be resolved without short baselines and radio flux 
would be missed, this in turn would increase the FIRC to 
higher values. In Fig.[B]we show the ratio of the flux-densities 
from NVSS and FIRST as a function of redshift. One can 
immediately see that at z < 0.2 the FIRST survey has miss- 
ing flux when compared to the NVSS. Such a selection effect 
could explain the recent results which sugge st an increase i n 
gm towards low (z < 0.5) redshift (e.g. iBourne et ail feoiO). 
Therefore, a combination of AGN contamination, and the 
resolving out of extended radio emission, could be respon- 
sible for the observed increase in giR towards low redshift 
based on deep field data which uses higher resolution radio 
imaging. 

However, if we return to the explanation in Section T3. II 
in which the longer wavelength e mission has an increas- 
ing contribution fro m evolved stars (I Rowan- Robinsonll2003l ; 
ICalzetti et aTfaldl , then we would expect the low-redshift 
sources to exhibit less correlation with the radio emission 
due to the fact that our data necessarily probe longer rest- 
frame wavelengths. This would not have a large effect within 
our dataset which is limited to z < 0.2 for direct detections, 
however when considering the evolution from z ~ 1 through 
to z ~ 0, this could explain the increase in giR towards the 
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Figure 8. The ratio of the flux-densities from the NVSS and 
FIRST survey as a function of redshift for radio sources in our 
sample with > 4cr detections in both surveys. Those sources iden- 
tified as AGN from their radio morphology are denoted by open 
stars. 

low redshifts, where there is an increased relative contribu- 
tion to the total far-infrared luminosity from dust heated by 
evolved stars. 

We therefore suggest that the evidence for an increase 
in the FI RC towards low red shifts found by previous au- 
thors (e.g. IBourne et al.l |2010| ) can be explained by one or 
a combination of resolving-out extended radio emission in 
high-resolution radio data coupled with an increased contri- 
bution from dust heated by evolved stars at progressively 
lower redshifts. 



3.2.2 Binning in far-infrared luminosity 

Figs.[6]and[7]show the giR and g2so parameters as a function 
of redshift, split into bins of far-infrared luminosity, respec- 
tively. We find no evidence that the mean giR or g2so in- 
crease with far-infrared luminosity, with all values between 
bins consistent within 1 a. However, it is also apparent that 
there is evidence for a negative trend in giR and g25o as a 
function of redshift in at least the most populated bin in 
Lfir at 10 < log 10 (ipiR/ L 0) < 11. In addition, we find 
that the stacked radio data provide values for giR ad g2so 
which closely resemble the values, within the uncertainties, 
on the individually detected sources, thus we can also be con- 
fident that we are not suffering any significant bias against 
sources falling below our 5(7 radio detection limit. 

One explanation for this negative trend with redshift 
is that any enhanced contribution to the radio emission 
from AGN activity in these sources would also decrease 
the value of gm and g250, and possibly increase the scat- 
ter. This obviously plays some part in the observed trend 
as the open circles in Fig. [6] denoting the sources with 
Lia GHz > 10 23 W Hz -1 , all lie below the mean values of gm 
and g250- For those sources with L1.4 ghz < 10 23 W Hz -1 
only a very minor enhancement of the radio flux would be 
needed to cause th e observed offs e t. Usi ng the radio contin- 
uum simulation of IWilman et all (|2008h we determine how 
many low-luminosity AGN would be expected over the red- 
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Figure 5. The FIRC as a function of redshift separated into three radio luminosity bins shown in the individual panels. The solid circles 
denote the 5a detections at the radio wavelength, the open stars denote those sources i dentified as AGN on the basis of jet-like radio 
structure and the open triangles denote those object selected at radio wavelengths from Hardcas tle et all ll201(]h . The solid horizontal 
line shows the weighted mean value of <jjr for all non-AGN sources with Li.4GHz < 10 23 W Hz -1 , and the dashed horizontal lines are 
the lcr standard deviation of the same sources. 



shift range < z < 0.2, where our sources with low values 
of Qir reside, and also with 10 20 < L1.4 ghz/W Hz -1 < 
10 . We find that we exp ect ^1 — 2 radio-quiet AGN or 
low-luminosity FR-I type (|Fanaroff fc Rilevl ll974T ) sources 
per square degree above the 5er radio flux-density limit 
of NVSS, thus in the SDP field we would expect 15-30 
sources where an AGN is contributing to the radio emis- 
sion at I/i.4GHz < 10 23 W Hz -1 . Thus it is entirely plau- 



sible that low-luminosity AGN could lower the value of 
giR for a small fraction of the sources in our sample. Fur- 
thermore, low-level AGN activity has been f ound to be 
relatively common in late- type galaxies (e.g. iFilho et alj 
2006; Pesroches fc Holl2009l ). Follow-up observations across 
several wavebands could be used to test this, e.g. high- 
resolution radio maps or deep X-ray imaging. 
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Figure 6. The FIRC as a function of redshift separated into three far-infrared luminosity bins shown in the individual panels. The solid 
circles denote the 5<r detections at the radio wavelength with Li.4GHz < 10 23 W Hz^ 1 and the open circles denote 5<r radio sources with 
Li .4GHz > 10 23 W Hz -1 . The open stars denote those sources identified as AGN on the basis of jet-like radio structure. We also show 
the results of our stacking analysis of the radio data at the position of the Herschel sources. The large filled squares and stars denote the 
average giR value from stacking at particular redshifts and far-infrared luminosities. In the top panel the filled square represents the mean 
9IR value at 8 < log 10 (LYiu/ Lq) < 9 and the filled star 9 < log lo (Lpiri/L0) < 10. In the middle panel the filled square denotes the value 
of qiB, from the stacked radio data for 10 < log 10 (LYiB,/ Lq) < 11 plotted at the mean redshift of the sources. Similarly in the bottom 
panel the filled square represents the stacked data at 11 < log 10 (Lpm/ Lq) < 12 and the filled star 12 < log 10 (LFiR/Z/Q) < 13 . The open 
squares and stars represent the same luminosity range but for radio stacks which give flux-densities such that Li.4GHz > 10 23 W Hz -1 , 
i.e. within the AGN-dominated radio regime. As in Fig. \E\ the solid line shows the weighted mean value of qir for each far-infrared 
luminosity bin, with the dashed lines showing the standard deviation of the data for Li.4GHz < 10 23 W Hz -1 . 
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Figure 7. <J250 as a function of redshift separated into three far-infrared luminosity bins shown in the individual panels. The symbols 
are the same as in Fig. [6] 



4 CONCLUSIONS 

We have used the H-ATLAS Science Demonstration Phase 
data to investigate the evolution of the FIRC over the 
redshift range < z < 0.5. Using a combination of the 
far-infrared data from H-ATLAS and the radio data from 
the FIRST and NVSS we find no evidence for evolution 
in the FIRC, with a median girt = 2.40 ± 0.12 over the 
redshift range < z < 0.5 . This is consistent with pre - 
vious work at z ~ (e.g. lYun. Reddv fc Condon! l200ll ). 
as well as studies focused on much higher redshifts (e.g. 



Ibar et al.l 120081; ISargent et all l20ld ; llvison et all l2010al lbl; 
Bourne et al.ll2010l ). in addition to empirical models based 



on radio and far-infrared luminosity functions (e.g. Wilman 
et al. 2010). 

Splitting the sample into bins of radio luminosity 
shows that in our far-infrared selected sample the FIRC 
is consistent with being constant with radio luminosity for 
sources restricted to the regime where star-forming sys- 
tems dominate the source counts, i.e. those with Li.4GHz < 
10 23 W Hz" 1 . We find that the dispersion in qm is a fac- 
tor of two higher than the dispersion for 5250- This could 
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be explained if the longer wavelength emission, i.e. in the 
sub-millimetre regime, was not produced from massive star 
formation and thus does not trace the same physical mecha- 
nism as the radio or the mid- to far-infrared emission. Such 
an effect has been found for lower redshif t galaxies where 
the SEDs have been studied in detail fe.g. lBendo et al1l20ld : 
ICalzetti et alll20ld) . This would imply that when calculat- 
ing star-formation rates using Herschel photometry it may 
be more accurate to restrict this to just using the shorter, 
rest-frame < 250 fim, data. However, we cannot entirely rule 
out the possibility of low-level AGN activity contributing to 
the radio emission and this would lower the value of qir. We 
also suggest that the increased dispersion in the value of qm 
could be linked to an increasing uncertainty in the models, 
which are less well constrained than the monochromatic flux 
at 250 fim, where we have data close to the rest-frame wave- 
length of interest. Deeper observations with PACS and/or 
observations at submm wavelengths could test which expla- 
nation is the more important. 

We also conclude that the increase in the ratio of far- 
infrared to radio luminosity towards low redshift found in 
previous work may be partly due to the relatively high spa- 
tial resolution of the radio data used. Such data would begin 
to resolve out extended emission from the large angular size 
galaxies which become more dominant in the low-redshift 
Universe. Future work combining the H- ATLAS survey with 
lower frequency Giant Metrewave Radio Telescope (GMRT) 
data at 325 MHz an d the Low Frequency Array (LOFAR; 
iMorganti et al.ll2010h survey data will also allow us to in- 
vestigate the FIRC over a larger range in redshift, due to 
the increase in depth of the radio data, and as a function of 
radio spectral index. This may provide important clues as to 
the link between the synchrotron emission traced by the ra- 
dio data and the thermal emission traced by the H- ATLAS 
data. 
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